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Abstract 

£Sj ' The loop expansion is applied to a chiral effective hadronic lagrangian; with the techniques of 

Infrared Regularization, it is possible to separate out the short-range contributions and to write 
them as local products of fields that are already present in our lagrangian. (The appropriate field 
variables must be re-defined at each order in loops.) The corresponding parameters implicitly 
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include short-range effects to all orders in the interaction, so these effects need not be calculated 
explicitly. The remaining (long-range) contributions that must be calculated are nonlocal and 
resemble those in conventional nuclear-structure calculations. Nonlinear isoscalar scalar (a) and 
vector (u) meson interactions are included, which incorporate many-nucleon forces and nucleon 
substructure. Calculations are carried out at the two-loop level to illustrate these techniques at 
finite nuclear densities and to verify that the coupling parameters remain natural when fitted to 
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the empirical properties of equilibrium nuclear matter. Contributions from the ujN tensor coupling 
are also discussed. 

PACS numbers: 24.10. Cn; 21.65. +f; 24.10.Jv; 12.39.Fe 



*Elcctronic address: \texttt{oberonjwm@yahoo . com} 

^Electronic address: \texttt{yihu_us@yahoo . com}; permanent address: Motorola Inc., 1501 W. Shure Dr., 
Mail Drop: 3-2G, Arlington Heights, IL 60004. 
■^Electronic address: \texttt{serot@indiana.edu} 



1 



I. INTRODUCTION 



In recent years, density functional theory (DFT) has emerged as a powerful tool for 
describing the structure of nuclei. Based on the seminal works of Kohn, Hohenberg, and 
Sham [1-3], DFT was originally developed for use in condensed matter physics; interestingly, 
it was discovered that it also provides a systematic framework for modeling the low-energy 
sector of the strong interaction [4-10]. Central to DFT is the notion that, if the energy 
functional is considered a functional of the density, then by minimizing the energy one 
can acquire the exact ground-state density [1, 3]. Furthermore, if the quantum many-body 
equations are replaced by single-particle equations with local, classical fields, one can exactly 
reproduce the total energy, scalar and vector densities, and chemical potential [2, 11]. Thus 
the many-body problem is reduced to finding the proper form of the energy functional. As 
it turns out, the lagrangian of quantum hadrodynamics determines just such a functional. 

Quantum hadrodynamics (QHD) is a low-energy effective theory of the strong interaction 
that uses the observed degrees of freedom at this energy scale (hadrons). It has evolved from 
the early model of Walecka [12] into more modern incarnations that are based on chiral effec- 
tive field theory (EFT) and DFT [7-10, 13-16]. The QHD energy functional is constructed as 
an expansion in powers of the mean meson fields, which are in fact Kohn-Sham potentials, 
over a heavy mass, such as the chiral symmetry breaking scale, and includes all possible 
terms consistent with the underlying symmetries of QCD. Field redefinitions are used to 
move the complexity of the problem (e.g., many-body forces and nucleon substructure) into 
the meson field self-interactions. As this is an effective theory, the terms are characterized 
by unknown coefficients. Once all the dimensional and combinatorial factors have been re- 
moved, however, the remaining dimensionless coefficients are assumed to be natural, or of 
order unity [7, 8, 17-20]. The long-range dynamics are included explicitly, and short-range 
physics is contained in the parametrization; therefore, this theory incorporates the natural 
separation of length scales. While QHD in principle contains all possible terms consistent 
with the underlying symmetries of the system, in practice it is a perturbative expansion for 
the energy functional that can be truncated at a manageable level (since the ratio of the 
mean fields to the heavy mass is small). The now finite number of coefficients are fixed by 
experimental data, and this theory can then be used for predictive purposes [21-27]. 

QHD is a strong-coupling theory. Unlike QED, there is no obvious asymptotic expansion 
to use to obtain results and refine them systematically. It is thus unknown whether QHD 
permits any expansion for systematic computation and refinement of theoretical results. 
One possibility is the loop expansion, which was partially explored in [28-31]. 

Why a loop expansion? The loop expansion is a simple and well-developed expansion 
scheme in powers of h that is derived from the path integral [28, 32-37]. The mean me- 
son fields are included non-perturbatively, and the correlations are included perturbatively. 
Therefore, one can analyze the many-body effects order by order. Previous work in QHD has 
shown that the mean-field terms dominate the nuclear energy, and exchange and correlation 
effects do not significantly modify the energy or nucleon self-energies, at least for states in 
the Fermi sea [34]. We stress that we are not certain that the loop expansion is practical; 
the answer to this question is left for future consideration. However, the loop expansion has 
the advantage that it is fairly easy to separate the short-range and long-range dynamics and 
to analyze their structures. 

The loop expansion was investigated in the simplest case at the two-loop level [28, 31] 
(only the most rudimentary couplings were retained). This was a severe truncation of the 
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EFT lagrangian and implies that we were interested only in the lowest level of accuracy. Here 
we are also concerned with the loop expansion only up to two-loop order (consideration 
of higher-loop effects is left for the future), but now nonlinear meson self-couplings and 
tensor terms are retained. Note that the nonlinear terms incorporate short-range many-body 
nuclear forces, vacuum dynamics, and the effects of nucleon substructure. The extension 
of the level of truncation to include these terms allows us to describe nuclear matter (and 
presumably nuclei) with greater accuracy. 

The result of the loop expansion up to the two-loop level is the mean field energy density 
[7-9] plus a series of integrals. We can categorize these integrals in two groups: baryon 
loops and purely mesonic loops. The baryon loops are integrals at the two-loop level that 
contain two factors of the baryon propagator and one factor of a meson propagator. Each 
nucleon propagator can be separated into a Feynman and a Density part [28, 34]. This is 
accomplished by incorporating the proper pole structure of the propagator. The Feynman 
part describes the propagation of a baryon or an antibaryon; the Density part involves only 
on-shell propagation in the Fermi sea and accounts for the exclusion principle. 

As before [31], the two-loop baryon integrals can each be separated into three distinct 
parts: the exchange, Lamb-shift, and vacuum-fluctuation contributions. The exchange 
term has two factors of the Density portion of the nucleon propagator. The Lamb-shift 
term contains both Feynman and Density parts of the nucleon propagator. The vacuum- 
fluctuation term contains two factors of the Feynman propagator. The Lamb-shift and 
vacuum-fluctuation parts are short-range physics, and it was shown in the simple case of 
[31] that they can be expressed as a sum of terms that are already present in the EFT la- 
grangian. This is an effective theory, so the coefficients of these terms are not yet determined; 
thus, the two-loop terms are just absorbed into terms already present in the lagrangian and 
are not to be calculated explicitly [31]. As a result, only the exchange portion of the two-loop 
integrals must be calculated explicitly. 

In addition, a number of purely mesonic loops arise when nonlinear meson field self- 
interactions are included. These terms, however, can be expressed as a power series in the 
meson fields with undetermined coefficients. As before, these terms are just absorbed into 
the lagrangian parameters and should not be calculated explicitly. 

It is interesting to note that the natural separation of the length scales in QHD and 
the parametrization of the short- distance dynamics is similar to Infrared Regularization 
[38-47]. In Infrared Regularization, loop integrals are separated into infrared regular and 
singular contributions (corresponding to short- and long-range physics). The regular portion 
is expanded as a power series in the momentum and is absorbed into the unknown coefficients 
of the underlying lagrangian; the singular part must be calculated explicitly. This is directly 
analogous to what was done for the loop expansion in QHD [31]. 

One question that remains unanswered is how the addition of these many-body correla- 
tions to the energy density affects the naturalness assumption, which states that once all 
the dimensional and combinatorial factors in a given term are accounted for, what remains 
is a dimensionless coefficient of order unity [7, 8, 17-20]. Once one decides on a level of 
truncation in the underlying lagrangian, the now finite number of unknown constants can 
be parametrized by the data. Thus, one purpose of this work is to prove that sets of pa- 
rameters exist at different levels of truncation in the two-loop energy density in which the 
parameters are all natural and thereby show that the naturalness assumption holds. For the 
purposes of this work, we will consider only uniform, symmetric nuclear matter, as in [31]. 
Here, however, we also include nonlinear meson self-interactions and tensor terms. These 
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additional terms not only increase the number of two-loop integrals, but also complicate 
the meson propagators, due to the short-range and many-body physics mentioned earlier. 
The loop expansion is carried out to the two-loop level, and the resulting two-loop integrals 
are separated into long-range and short-range physics. The short-range contributions are 
expressed in forms that already appear in the lagrangian. Since the coefficients of these 
terms have yet to be determined, they are just redefined, thus incorporating these new con- 
tributions into the mean field energy (as in [31]). As a result, they are already present in 
the one-loop QHD calculation. The long-range physics must be calculated explicitly. In this 
work, we fit the two-loop energy to empirical properties of nuclear matter. Then we compare 
parameter sets constructed at the two-loop level with those of previous works obtained at 
the mean field level and consider the naturalness of the parameters. 



II. THEORY 



A. Isoscalar Meson Self-Interactions and Mixing 

In this section, we investigate the effects of nonlinear meson self-interactions (including 
isoscalar meson mixing) on the loop expansion. Here we will keep terms up to order v = 4 in 
the isoscalar meson sector of the QHD lagrangian [7, 8] (i.e., k 3 , k 4 , rji, r] 2 , and Co 7^ 0), and 
let ai = a 2 = 0. Since we will calibrate the parameters using nuclear matter properties only 
(and not finite nuclei), we want to keep the number of parameters at a manageable level. 
Therefore, the baryon and pion parts of the lagrangian are unchanged from [31] (except for 
the addition of the two-pion-nucleon vertex): 



C = -if; 
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7m {dp - W - i-rlnlbd^TL + — 7 M [vr, <9 M vr] + (M — g s <f>) 



-\ ( W - \« , (i) 

which follows from the chirally invariant lagrangian in [7] by retaining only the relevant 
(lowest-order) terms in the pion fields. 1 The neutral scalar and vector meson lagrangian is 
now 

r 1 fa a\ 2 2,2 f 1 , K z9s<t> .K A g 2 s 4> 2 \ 1 

C *» = -2 {d ^ ~ ms<P \2 + 3\W + 4! M^J~4^^ 

-\mlV^ (l + ^ + | ^) + l - C o9l ( V^f , (2) 

where V^ u = d i y v — d v V iJi and n = \rc a -T a . Here tp are the fermion fields and 0, V^, and 7i a are 
the meson fields (isoscalar-scalar, isoscalar- vector, and isovector-pseudoscalar, respectively). 
The heavy meson fields are also chiral scalars. Note that in this work, the conventions of 
[10] are used. As in [31], we are not making a chiral expansion in powers of the pion mass 
and include it for kinematical purposes only. 



In Appendix A, we discuss the role of the tensor coupling to the neutral vector meson at the two-loop 
level. 
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The generating functional is defined in the usual way by 
Z\j,J„] = exp{tW[j,J,]/h} 

= AT 1 [ D^D^D^DiV^Di-Ka) 



X 



cxp |l j d 4 x [C(x)+j(x)<j>(x) + J,(x)V,(x)]} , (3) 



where J\f is the normalization factor (in effect, the vacuum subtraction, which is also an 
expansion in loops), j(x) and J^{x) are the external sources corresponding to the meson 
fields and V^, respectively, and the connected generating functional is W[j, J M ]. The path 
integrals of interest here are those over the isoscalar quantum fluctuation fields a and 77^ 
(the others are discussed in [31]). With %, r] 2 ^ 0, the integral is bilinear in form, or 



1 = J D(a)D(rjr) exp ji J d 



x 



-aA s 1 a+ -fj^D *fj v 
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M >' M- 

where the meson propagators in momentum space are 
A-i (k) _ P 2 / gsfo <u gUl \ _ 2 HiAv v° 

A 5 (fcj - k + m s ^L + k 3 M + 2 M2 J m y 2 ^ ^ ^ , 



(4) 



(5) 



2 M 2 
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A - k k _ T/°T/° 



(6) 



To simplify the equations, we introduce a five-component notation for the meson fields 
and their corresponding auxiliary sources [48]: 



s = 



a 



u 

U a 



(7) 



where s a and t a have a — 0, 1, . . . , 4. We can then rewrite the path integral as 



J D(s a ) exp J d 4 
= exp | -- tr In [A" 1 ] 1 exp 1 t a A al3 tp 1 , 



where the matrix A is 



and 



A 5 1 (A;) V V U ° 
VV° V~l{k) 



2 / Vl9S . , \ 



(8) 



(9) 



(10) 
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If rji — 7)i — 0, then A' 1 is diagonal, there is no mixing, and the system reduces to the 
non-mixing case, a simple example of which was discussed in [31]. 

Now we want to solve for the components of A(k) using A(k)A~ 1 (k) = 1: 



^00 A 0n 

A u 'o A v i^ 



vv° v-l{k) 

This leads to the following four equations: 

A m [^- S \k)]+A^ = 1, 
A m r,V? + A,,[V-l{k)\ = 0, 
A V , Q [AgWj+A^vVZ = 0, 
A v , riV° + A u/fl [D-i(k)] = <W . 

Solving these, we get (V 2 = V°V°) 
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k 2 + M 2 V - %g 2 v V 2 
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k 2 + Ml- ^g 2 v V 2 



x 
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M'y ~ ^ 9 yV 2 
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C, 



Co 



—r)[k z + M' v — jg 2 v V 2 



4,0 = -?? ( A; 2 + M 2 - |<£y* 



M 2 - |^V 2 ) V° + (k -V) k u , 



(11) 

(12) 
(13) 
(14) 
(15) 



(16) 



/ C , (17) 
/C , (18) 



Av'fX 



where 



k 2 + My — ^g 2 v V 2 \ (m 2 v - ^g 2 v V 2 \ - C jg 2 v (k ■ Vf 



(k 2 + M|) 

+ (k 2 + M 2 V - C jg 2 v V 2 ^ k»k v , + [m 2 - C jg 2 v V 2 ^ VX> 
+ ( {g 2 v (k-v) {kX' + k^) 



M 2 V - C jg 2 v V 2 )V 2 + (k- V) 2 



<W + V 2 k il k v , 



- (m 2 - C fg 2 v V 2 ^ VX -(k-V) (kX> + k.Vl) }) /C , 



(19) 



C = 



k 2 + M 2 y~ ^ 9 2 yV 2 
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x \(k 2 + Ml) 



M 



v 



k 2 + M 2 v - C jg 2 v V 2 



^V 2 ) V 2 + (k ■ Vf 



Ml 



Co 2/7 T/\ 2 



and 



Ml 



Ml 



= m s [l + K 3 



9s<Po , lugjn 
2 M 2 



M 



+ 



m 2 ^f^ 2 



= m y 1 + 771 



M 2 M 2 



(20) 

(21) 
(22) 



The corresponding path integral in the vacuum subtraction has a similar structure, except 
with vanishing classical meson fields; therefore, A is replaced by A , or 



A (k) = 



A° s (k) 

vim 



(23) 



where and TP are the free propagators for the scalar and neutral vector mesons, respec- 
tively: 



A°,(A:) = 



1 



k 2 + m% — ie 



k 2 + my — ie 



k k 



m 



V 



(24) 
(25) 



The pion self-couplings do not contribute to the energy density if there is no pion con- 
densate. Therefore, we simply ignore the pure pion terms in the following calculation. The 
remaining generating functional is 



Z\j,J li ]=tf- 1 exp\j: J d 4 x [C (x) + j{x)<j>o{x) + Ux)V°(x)} 



x exp 



tr In [G° F (k)G- H \k)} - - tr In [A (k)A~ 1 (k)] 



x 



exp \ iH 1/2 / d 4 x 



, ■ 9a a 



iS 



9s 



-iS 



8£(x) J \ L<5m(x 
-iS 



-iS 



' -iS ' 






■>)[ 



-iS 



+ my 



Vi9s V29s , 
2M 2M 2<Po 

-iS 



5U^x) 



-iS 
8u(x) 

-iS ' 

5U u (x)_ 



5u(x 
-iS 



5U,(x) 

-iS ' 

SU v {x)_ 



1 3 

-iS 
WJx) 

2 mg 2 s v o 



+ m 



v 



2M 2 M 



-iS 
Su(x) 



1 2 



-iS 

mjx) 



7 



+ h / d X — ( g v 



24 



-iS 



m 



2 »Ms 



V 



4M 2 



x exp < i 



5u(x) 
i J d 4 x J d 4 y 



-iS 
[5Ujx) 

1 

2 



-iS 



-iS 



5U v {x) 



-iS 



5U u (x) 



— m 



2 K ^9s 
S 2AM 2 



-iS 



1 4 



Su(x) 



5U^x)_ 
t a A a/3 (x - y)tp 



+ - Ca(x)Af(x - y)( b (y) - £(x)G H (x - y)£{y) 



(26) 



-J sources = 



Here ( a is the pion source and £ is the nucleon source. We define the contribution at the 
mean field level as 



1 



C = -^mlvX[l+ Vl ^ + 



2 " "v ■ fj, ji 



9S<Po , ^2 9sK 



M 2 \P ' 



+ 



+ 



2 3! M 4! M 2 



and the pion and fermion propagators are given by 

1 



G° F (k) 
G H (k) 



k 2 + m 2 — ie 



— Sab , 



-1 



ift+M — ie ' 
-1 



(27) 

(28) 

(29) 
(30) 



It is straightforward to acquire the connected generating functional at the one-loop order 

d 4 k 



and therefore 



rW [0 e , K°] = / ^ {A) + ih J ^ (tr In [(^(AOGhW] 



- - tr ln [A (k)A-\k)] 



(32) 



The baryon term is identical to the simple case in [31]. To treat the isoscalar terms, we first 
consider the identity 



tr ln [A (k)A~ 1 (k)] = lndet [A (k) A~ 4 (k)] , 
since it is easy to acquire the determinants: 

det [A (k)} = \m\ (k 2 + ml) (k 2 + m 2 v ) 3] 



(33) 



(34) 
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det [A~\k)} = [e + M v - C jg v V ( 



x 



(m v - |^k 2 ) + f [k 2 v e 2 - (k ■ v e f] 



k 2 + M v - C jg v vf 



-v 



M 2 V -%g 2 v vnvr (k-V 



(35) 



where M| and M v are denned by Eqs. (21) and (22), and V° {1) = V^ 6 = i5^V 
Apparently, det [A _1 (/c)] is an analytic function of e and V e , and 

det [A-\k)] (0 e = 0, V e = 0) = det [A„ X (A;)] . 



(36) 



Therefore, the meson contribution to the effective action can be categorized by the general 
form 



(37) 



i,j = 



where the prime indicates that % + j > 0. However, terms of this type already appear in 
the general QHD lagrangian; since the coefficients of these terms have yet to be determined, 
these meson contributions can be absorbed into the preexisting terms. As a result, they 
don't need to be separately considered. 

Now we want to determine the connected generating functional at the two-loop level, 
For this, we have to make a Legendre transformation to get the effective action. The 
classical meson fields are approximated as [31] 



/ <SWi 



<j> (x) 



</> e (x) 



+ / d 4 yA(x - y) 



5 My) 

5W l 



\ 



J 



+ 0(h 2 ) . 



(38) 



One can verify that the tadpole diagrams in W 2 eventually cancel, and so we are left with 
only the one-particle irreducible diagrams: 



r (2) [0 e ,v; e ] = 

— ih f d 4 x 

+ h 2 I d'-.v 



d 4 xC [Ux),V^x)} 

d 4 k 



tr In [Gp(k)Gff 1 (k)] - - tr In [A (k)A~ 1 (k)] 

2 



(2^) 4 
d 4 k 

(2tt) 4 J (2tt) 4 [ 2 



d A q / gis 



: -A 00 (k-q)tr[G H (k)G H (q)} 



9 2 v 



- ^A^k - q) tr [G H {k) llx G H {q) lv \ 



+ i9s9vA 0ll (k - q) tr [7 ^G H (k)G H (q)] 



j±A*(k-q)tr [G H (k)(V- <t)l^G H {q){U-- d) l5 j 
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> , 9sm 2 s g 2 s m 2 s (j) 



^^-A 00 (k)A 00 (q) 



A m (k)A 00 (q)A 00 (k - q) 



+ -r] 2 [A 00 {k - q)A^(k)A lxu (q) + 2A 0fl {k)A 0u {q)A^{k - q)} 



9sm 2 s 
M 



+ K A 



9s m s^ 
M 2 



A 00 (k - q)A 0ll {k)A Qtx {q) 



\ A M A M + 2A Q ,(k)A Q ,(q)} 



+ l A ^( k ) A M + 2VC=)^(?)] 



+ 



2M 4 



v 



[A 00 (k)A 00 (q)V*A^(k - q)V: + 2A 00 (k)V^A IM0 (q)A 0u (k - q)V:] 



m 



2 V2(o9s9y_ 
6M 2 



v 



[A, x {k)A x ^q)V^A u0 {k -q) + 2A 0X (k)A X0 (q)V*A llu (k - q)V v e ] 



1 



+ - Co9v [A, u (k)A, v (q)VZA xx ,(k - q)V^ 

+ 2V;A^k)A vX {k - q)A xx ,(q)VZ] - VEV \ , 



(39) 



where VEV is the vacuum subtraction. Notice that there are now four distinct fermion loops 
at the two-loop level. These reduce to the familiar three loops from the non-mixing case 
when i]i = r] 2 = 0. 

With dimensional regularization, we can shift momentum variables to remove the depen- 
dence on the vector field from the fermion propagators and thereby replace Gu(k) by G*(k), 
where 



G*{k) 



-1 



i^ + M -gs<po 



_k 2 + M* 2 -ie E*(k) 
= G* F (k) + G* D (k) . 



S[h- E*(k)]9(k F -\k\) 



(40) 



Here M* = M — g,s<Po- We can then separate each of the fermion loops [using G*(k) = 
G* F (k) + G* D {k)} into three contributions: an exchange term, a Lamb-shift term, and a 
vacuum-fluctuation term. The exchange terms will be evaluated explicitly later. For now, 
let us consider the Lamb-shift term for the fermion loop involving scalar-vector mixing: 



^ A ^-^b,G* F {k)G* D {q)\ 
f d 4 k f 



d(1 77T 



M*(k + q), 
k 2 + M* 2 
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x [(k - qf + M v 2 ] {[(k — q) -V](k — q), + M* 2 V») 
[(fc _ g ) 2 + M 2 ] [ {k - qf + My 2 ] 



x 



x 



(* - g) 2 + M* 2 + ^g 2 v V 2 



M, 



v 2 -^[(*-?)-V] a 



oc 



- V 2 ([(k ~ qf + M* 2 ] M* V 2 + 
d A q f d 4 p 6 (kf 



q 



(* - g) 2 + M* 2 + |^\/ 2 
6(qt-E*(q)) 



[(k-q)-Vf 



(2tt) 4 7 (2tt) 4 £*(g) [p2 + 2 (p. g ) + g 2 + M *2j 
x [p 2 + My 2 ] [{p .V)(p 2 + 2p-q + M* v 2 ) + 2M* 2 (g • V)] 



x 



P 2 + My 2 + ^g'yV 2 



M, 



-v 



( P 2 + Mj) ( P 2 + My') 

( P 2 + My') My'V 2 + (p 2 + My 2 + J 9^ ( P • V) 



2 Co 2 f t/\2 

- t9v (p ■ V) 



-1 



oc 



i = o i = o 

21 i oi«2 



x(p.g) l (p.^{(p.\/)[p 2 + 2(p-g) + M* 2 ] +2M* 2 (g-\/)} 



oc 



y 



d 4 g (Jfe F - |q| 



) 4 £*(g) 



5(g 4 -£*(g))F(M*,\/ 2 ,g-\/) , 



where 



2 _ ^2 to^yl 



My^ = My — 
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(41) 



(42) 



kp is the Fermi momentum, and E*(q) = (q 2 + M* 2 ) 1 / 2 . We have made the substitution 
Pn — (k — q)n-> which is allowed if we keep q fixed in the k integral. The function F comes 
from the fractional integrand and is an analytic function of the variables p 2 , g 2 , V 2 , and M*. 
The function F comes from the integral over p and is analytic in the variables M*, V 2 , and 
q ■ V, since p^ is always spacelike. The on-shell condition g 4 = E*(q) is used to eliminate 
the dependence on g 2 in the last line. 

Remember that V® = iVoS^. Thus, after integrating over g 4 , the remaining expression 
is an analytic function of E*(q). The final result will be complicated, but should be explic- 
itly dependent on some types of baryon densities. These terms can be directly related to 
terms in the effective lagrangian that are some combination of the nucleon fields (with some 
derivatives) multiplied by some Lorentz-invariant polynomials in the mean meson fields; as 
a result, they can be parametrized by local terms in the effective lagrangian. A simple case 
was illustrated in [31]. 

The term with the pion propagator is the same as in the non-mixing case. The fermion 
loops with A 00 and will, in the case of r\\ = r] 2 = 0, reduce to the equivalent loops 
with the pure neutral scalar and vector mesons from the non-mixing case. We expect that 
they will behave similarly in both the mixing and non-mixing cases, with the former simply 
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leading to more general local combinations of the fields. Thus the Lamb-shift terms can be 
parametrized by terms already present in the lagrangian (before truncation) . 

The vacuum-fluctuation terms and all the pure meson loops have no explicit density 
dependence, so all the conclusions of the previous analysis [31] are still applicable here. 
The short-range physics which is contained in the Lamb-shift and vacuum-fluctuation terms 
is absorbed in preexisting terms in the QHD lagrangian. The long-range physics arising 
from the exchange term must be explicitly calculated; these integrals are considered in the 
following section. 



B. Two-loop Integral Formulas 

We first consider the integral with the propagator from Eq. (16). Therefore, one can write 



where 



(2) 

4>-EX 



19s 



1 d 4 k r d\ 

kF |k| 2 rf|k| r kF \q\ 2 d\q\ 
E*(k) J 



32^ ./ E*(k) J E*(q) j_, 



d{cos 6) 



x [E*(k)E*(q) - |k||q| cos# + M* 2 ] 

x \2E*(k)E*(q) - 2|k||q| cos# - 2M* 2 + My 2 ] 



x 



2E*(k)E*(q) - 2|k||q| cos# - 2M* 2 + My 2 + —g v V 2 



M, 



v 



B 



-^g v V 2 [E*(k)-E*(q)] 2 y B , 

\2E*(k)E*(q) - 2|k||q| cos# - 2M* 2 + M 2 ] 
x \2E*(k)E*(q) - 2|k||q| cos# - 2M* 2 + M* 2 ] 



x 



2E*(k)E*(q) — 2|k||q| cos# — 2M* 2 + M v 2 + —gyV 2 



C °-g 2 v V 2 [E*(k)-E*(q)} 2 



M 



* 2 
V 



(43) 



+ r] 2 V 2 {[2E\k)E\q) - 2|k| |q| cos # - 2M* 2 + M£ 2 ] M; 
2E*(k)E*(q) - 2|k||q| cos9 - 2M* 2 + M* 2 + — g v V 2 



x[E*(k)-E*(q)] 2 ) . 

Now we consider the integral with the propagator from Eq. (19); as a result, we get 

9± f d*k 
2 J (2tt)4 J (2tt) 



(44) 



?(2) 

'V-EX 



%A„ u (k-q)trty,G* D (k)l„G* D (q)} 
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Wv ["' IM 2 rf|k| {»' |q| 2 d|q| /•' 

x [[2£*(fc)£*(g) - 2|k||q| cosfl - 2M* 2 + M|] 



x <^ 2 



2E*(k)E*(q) - 2|k||q| cos# - 2AT 2 + My 2 + ^-g v V£ 

- C jg 2 v V 2 [E*(k)-E*(q)] 2 
x - |k||q| costf - 2M* 2 ] 

+ ^g 2 v V 2 M v 2 [E*(k)E*(q) + |k||q| costf + M* 2 ] | 

-rfV 2 (2{[E*(k) - E*{q)f - M* 2 } [E*(k)E*(q) - |k||q| cos# - 2M* 2 ] 
—My 2 [E*(k)E*(q) + |k||q|cos# + M* 2 ])]/5 . (45) 

Lastly, we consider the integral with the mixed propagator from Eq. (17), which becomes 

e m-ex = -igsgv / j -^^^^-^^h^D^G^q)} 

jg s gvV M *Vo 



M y 2 



2 f k * \k\ 2 d\k\ [ k * \q\ 2 d\q\ f 1 



l(i;r ; 

x [E*(k) + E*(q)} [2E*(k)E*(q) — 2|k||q| cos# — 2M* 2 + My 2 ] j B . (46) 
The formula for the pion contribution to the two-loop energy density is the same as in [31]. 

III. DISCUSSION 

In the previous section, we showed that the short- and long-range dynamics of the two- 
loop contributions can be separated even in the presence of various nonlinearities in the 
isoscalar meson fields. The short-range physics can be expressed as a series of terms that 
are already present in the underlying lagrangian; thus, they are implicitly contained in 
the undetermined coefficients. The long-range exchange contributions must be calculated 
explicitly. The question remains what effect these new contributions have on the naturalness 
assumption. 

Naturalness [7, 8, 17-20] states that once all the dimensional and combinatorial factors in 
a given term of the lagrangian are removed, what is left is a dimensionless constant of order 
unity. To show that naturalness is not violated by the inclusion of the two-loop exchange 
integrals, we will show that sets of natural parameters exist that adequately reproduce the 
empirical properties of symmetric nuclear matter at equilibrium [13]. 

We conduct parameter fits at different levels of truncation in the underlying lagrangian, 
as in [7], but here we include the two- loop exchange contributions. Note that the two- 
loop portion of the energy density contains no new parameters. In this work, we fit the 
parameters to bulk properties of nuclear matter: £ / p B — M = —16.10 MeV, kF = 1.30 
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1.7424 


1.2072 


1.6698 


4.3207 


2.6308 


2.1978 


K4 


-8.4836 


-1.8824 






-5.6463 


-5.3909 


Co 


-1.7750 


9.9968 






8.9212 


9.0313 



TABLE I: Parameter sets used in this work compared with one-loop parameter sets at the same 
level of truncation. 



fin -1 , M* /M = 0.610, and compressibility K = 250 ± 50 MeV. A least-squares fit was 
conducted to these four pieces of data with respective weights of 0.0015, 0.002, 0.0015, and 
0.08. The fitting was performed using the downhill simplex method, with nuclear equilibrium 
maintained throughout the procedure. Table I lists the results of this fitting for a number 
of different levels of truncation in the lagrangian (MIA, M2A, M3A, and M4A) — defined 
below — along with some corresponding parameter sets determined at the one-loop level with 
the same truncation (Q2 and CI [7]) for comparison. We hold the meson masses fixed in the 
new sets (ms/M = 0.5400 and my/M = 0.8328), since there are only four empirical input 
data; they accurately reproduce the reduced meson mass in [49]. We let M = 939 MeV. 

• The mean field lagrangian for the MIA (and Q2) set includes cubic and quartic scalar 
field self-couplings, as well as the quartic vector field self-coupling. In this case there is 
no mixing in the meson propagators (77 = 0); as a result, there are only three two-loop 
exchange integrals (A 0fl drops out). 

• The M2A (and CI) set includes a cubic scalar field self-coupling in addition to the 
third-order mixing term. This case contains mixing in the meson propagators through 
the inclusion of the rji term; therefore, there are four exchange integrals to calculate. 

• Two other mixing sets were constructed (M3A and M4A); however, no comparable 
mean field set exists in the literature. These sets have the same terms as MIA plus 
r/i 7^ and f| 2 ^ 0, respectively. Both sets give rise to four two-loop integrals as in 
the M2A case. 

In the set MIA, the value of Co was allowed to vary; however, Co places a strong constraint 
on the compressibility K. One must also ensure that the vector field equation still has an 
appropriate solution at high density. Confining Co to positive values is sufficient to reproduce 
the correct high-density behavior. One can see from Table I that parameter sets exist at 
the two-loop level in which all the parameters are natural. In all cases the properties of 
nuclear matter were reproduced accurately while maintaining naturalness. The bulk nuclear 
properties for the various sets are shown in Table II; note that all the observables in this 
table except gyV were used in the fit. 
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0.532 


0.611 
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292 
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255 
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K 


279 
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235 
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TABLE II: Nuclear matter properties for the parameter sets in Table I. The values of £/pB ~ M, 
<7yVo, and K are in MeV. The Fermi wavenumber is in fm -1 and M*/M is dimensionless. 





Q2 [7] 


MIA 


CI [7] 


M2A 


M3A 


M4A 


c(2) 

°4>-EX 


38.77 


33.51 


26.73 


21.26 


31.39 


32.99 


c(2) 

°V-EX 


-28.32 


-15.49 


-23.47 


-18.50 


-15.85 


-15.67 


c(2) 

°M-EX 






3.54 


11.99 


4.24 


3.30 


c(2) 


12.87 


12.76 


12.73 


12.79 


12.71 


12.68 



TABLE III: Size of two-loop integrals for the parameter sets in Table I. Values are in MeV. 



The size of the various two-loop integrals is shown in Table III. We also show the size of 
these integrals when the constants are taken from the Q2 and CI sets. Note that there is 
some cancellation between the scalar and vector contributions. 

Next, we compare the magnitude of the two-loop integrals to the mean-field terms. In 
Fig. 1, the mean field terms for the Q2 and MIA sets are plotted order by order in the 
power counting, with the two- loop integrals calculated using MIA. The crosses in Fig. 1 
represent the expected magnitude per order of terms in the meson sector using the rules of 
naive dimensional analysis (NDA) [7], with the chiral symmetry breaking scale of A = 650 
MeV. Fig. 4 compares the saturation curve of MIA at the two- loop level with the saturation 
curve of Q2 at the one-loop level and the two-loop binding curve of Q2. Apparently, the net 
effect of the exchange terms is repulsive. 

In Fig. 2, the mean field terms for the CI and M2A sets are plotted order by order in 
the power counting, with the two-loop integrals calculated using M2A. A comparison of 
the saturation curve of M2A at the two- loop level with the saturation curve of CI at the 
one-loop level and the two-loop binding curve of CI is shown in Fig. 5. 

Lastly, the terms (mean field and two-loop) for both M3A and M4A are plotted in Fig. 
3. Fig. 6 is a graph of the saturation curves for both M3A and M4A at the two- loop level. 

One can see from Figs. 1, 2, and 3 that the magnitudes of the mean field terms are 
consistent with the NDA estimates to within a factor of two, which supports the naturalness 
of the coupling parameters. 2 Moreover, the magnitude of the two-loop terms is shown 



2 The determination of naturalness is a tricky business, and the relevant observation here is that the 
hierarchy of successive terms in Figs. 1, 2, and 3 is maintained with the new parameter sets. A more precise 
discussion of the naturalness would require an examination of the linearly independent combinations of 
parameters, as in [13]. 
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FIG. 1: Comparison of the magnitudes of the mean field terms in the meson sector with the two- 
loop exchange integrals. The inverted triangles represent, from top to bottom, the scalar, vector, 
and pion two-loop integrals. The abscissa denotes the order v in the power counting at the mean 
field level [7]. 



to be roughly equal to the third order [y = 3) in the power counting in the mean field 
energy. While this is not large, they cannot be neglected in a description of nuclear matter 
properties, particularly in view of the nearly complete cancellation of scalar and vector terms 
at order v — 2. Nevertheless, when the parameters are properly renormalized to reproduce 
empirical equilibrium properties, exchange corrections modify the energy only slightly, and 
the two-loop binding curves for nuclear matter are similar to the one-loop results. 

Exchange interactions are two-body, so they go like 0(p 2 ) in £ . This is the same depen- 
dence as the two-body mean field terms at 0{y — 2). Because of the lack of a spin-isospin 
sum in the exchange terms, the exchange terms are numerically smaller and contribute at 
the same level as the 0(u = 3) terms. A deeper understanding of how loops fit into the 
finite-density power counting scheme will require a look at higher-order terms in the loop 
expansion, which is currently under investigation [50]. 

Exchange contributions from the the vector meson two-loop diagrams with tensor cou- 
plings at one or both vertices were also considered. The analytic expressions for these 
integrals are given in the Appendix, and the numerical results are shown in Table IV. The 
vector-tensor cross terms (XT) are roughly equal to fourth order [y = 4) in the power count- 
ing, and the tensor-tensor terms (TT) are negligible. This is not surprising considering that 
the tensor coupling on a vertex introduces an additional factor of V/M ^ 1/4 (see [7] for a 
description of the power counting). A value of fy = 0.5 was used here, which is comparable 
to values obtained from fits to properties of finite nuclei [7]. The tensor terms are small 
enough that their effect on the naturalness of the parameters is negligible. 

In conclusion, we have shown that the results of [31] are still valid when various nonlinear- 
ities in the isoscalar meson fields are included in the lagrangian. The short- and long-range 
physics at the two-loop level in the loop expansion can be separated. The short-range (local) 
dynamics is parametrized by the undetermined coefficients of the lagrangian, like those in 
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FIG. 2: Comparison of the magnitudes of the mean field terms in the meson sector with the two- 
loop exchange integrals. The inverted triangles represent, from top to bottom, the scalar, vector, 
pion, and mixed two-loop integrals. The abscissa denotes the order v in the power counting at the 
mean field level [7]. 





Q2 [7] 


MIA 


Cl [7] 


M2A 


c(2) 

°TT-EX 


-0.0241 


-0.0189 


-0.0210 


-0.0195 


c(2) 

C XT-EX 


-3.94 


-2.95 


-3.43 


-3.12 



TABLE IV: Size of two-loop integrals involving the tensor coupling to the vector meson for the 
parameter sets in Table I. Values are in MeV. 



the nonlinear meson self-interactions, and is either removed by field redefinitions or fitted 
to empirical data. The long-range (nonlocal) physics must be calculated explicitly. The as- 
sumption of naturalness is shown to hold in the presence of the two-loop exchange integrals. 
The size of the two-loop contributions is roughly equal to third order in the mean field power 
counting. The inclusion of tensor contributions was also considered and found to be small. 
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FIG. 6: Comparison of the saturation curves for the sets M3A and M4A (both at the two-loop 
level) . 



APPENDIX A: TENSOR PORTION OF TWO-LOOP CALCULATION 



In this section, we present the corrections to the theory arising from the tensor coupling 
to the neutral vector meson at the two-loop level, where 
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After working out the variational derivatives, we acquire the energy density. Now we can 
replace the Hartree propagator by Gh — > G* = G* F + G* D , as before. The Lamb-shift and 
vacuum fluctuation terms are absorbed by the same process as before. (The only difference 
here is more gamma matrices and they don't influence the result.) Now we are left with 
only the exchange terms 
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which were reduced in the same manner as in the previous discussion. Here B is defined in 
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